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O . Abstract 
O 

, This paper establishes and tests procedures which can determine the electron 

^ ' energy gap of the high-temperature superconductors using the t — J model 

, with spinon and holon quasiparticles obeying fractional statistics. A simpler 

Q , problem with similar physics, the spin susceptibility spectrum of the spin 1/2 

O 

fractional-statistics gas, is studied. Interactions with the density oscillations 
of the system substantially decrease the spin gap to a value of (0.2 it 0.2) 
huJc, much less than the mean-field value of hcoc- The lower few Landau 
levels remain visible, though broadened and shifted, in the spin susceptibility. 
As a check of the methods, the single-particle Green's function of the non- 
interacting Bose gas viewed in the fermionic representation, as computed by 
the same approximation scheme, agrees well with the exact results. The same 
mechanism would reduce the gap of the t — J model without eliminating it. 
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I. INTRODUCTION 



In this paper, we investigate the excitation spectrum of the spin ^ fractional-statistics 



gas, focusing upon the minimum energy required to flip a spin. We develop techniques 
enabling the computation of the energy gap of the cuprate superconductors in terms of the 
t—J model possessing spinon and holon quasiparticles obeying fractional statistics Jll . Using 
that framework, Tikofsky, Laughlin, and Zou determined the optical conductivity and 
other quantities in good agreement with exact diagonalizations of the t — J Hamiltonian for 
small systems. However, their computed electron energy gap is roughly five times larger 
than the experimental limit 0] for the cuprates. In addition, their linear response functions 
displayed unreasonable sharp structure, which had to be smoothed out in an arbitrary 
fashion. The spin-one excitation spectrum of the fractional-statistics gas provides insight 
into the electron gap of the more complicated t — J model. Specifically, an electron hole 
in the t — J model consists of a composite of two particles, a spinon and a holon, each 
obeying fractional statistics, just as the spin-one excitation of the fractional-statistics gas 
is a composite of a spin-up particle and a spin-down hole. As will be demonstrated later, 
interactions with phonons, the density oscillations of the fractional-statistics gas, broaden 
the spin susceptibility spectrum, thereby lowering the spin gap, the minimum spin-one 
excitation energy. As a possible comparison, Girvin and collaborators found the exact 
ground state of the spin | fractional-statistics Hamiltonian with an attractive delta-function 
interaction. However, they neglected the broadening of the spin susceptibility spectrum and 
consequently obtained an overly large spin gap. 

Let us define the system being studied. The non-interacting fractional statistics gas Q 
with spin corresponds to the spin S fermionic Hamiltonian, 



anyon 





where the gauge field Aj has the value 
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and Fj denotes the two-dimensional vector locating the ith particle. In the problem being 
studied, the particles have spin S = ^ and statistics = I- The cutoff of the gauge field A 
at small distances prevents a weak spin instability investigated by Beran and Laughlin 0. 
Setting the cutoff parameter a equal to 10'^ avoids the Hartree-Fock instability described in 
that paper with a negligible effect upon the subsequent spin gap calculations. In the lattice 
system to which the fractional-statistics formalism will be applied, the gauge force would be 
cut off at a distance of one lattice spacing. 

The spin susceptibility, defined for an unpolarized spin | system to be 

X(ri h I r2 ta) = < $ I T{a\r^ t^) a\r^ ta)} I $ > 

= -2« < $ I T{(T- (ri ti) a+(r2 h)}\<^> , (1.3) 

determines the spin gap of the fractional-statistics gas. The wavefunction |$> is the ground 
state, and the operator time dependence follows the Heisenberg picture, 

d{t) = e^i/mtQ^-m)Ht _ ^^ 4) 

The 2;-component of spin is 

a^(r) = ^|(r)x[/^(r) - ^I(r)xl>;(r), (1.5) 
while the spin-raising and spin-lowering operators are defined by 

a+{r) = ^\{r)^^{r) (1.6) 

and 

a-(r) = vl/t(r)vl/^(r), (1.7) 

where ^J.(r) and ^E'o-(r) are fermionic operators creating or destroying a spin a particle at 
the position r. The Fourier transform of the susceptibility, 

X(Q, ^)=L'JJ e-^(Q-^--*) x(r t \ 0) dt dr , (1.8) 
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determines the spin gap A, the smallest positive frequency uo for which Im ^(Q, ^) is non-zero 
at some momentum Q. 

Finding the spin gap at the mean-field level is straightforward. The system may be 
approximated by the mean-field Hartree Hamiltonian Tio, 



1=1 



^ 1 ^2 



Pi + -A, 

c 



:i.9) 



corresponding to particles in a constant magnetic field 



where 



A = ^Bxr, (1.10) 



he 

B=(l-z/)— pz (1.11) 

e 



and p is the average particle density of all spins. The resulting energy spectrum consists of 
sharp Landau levels, with a gap between Landau levels of the cyclotron energy, 

heB , , ^ 

huc = = 27r(l-z/)— p , (1.12 

mc m 

as displayed in Fig. |I[ For spin | semions, the lowest Landau level is occupied with particles 
of both spins, and the other levels are empty. Thus, the spin-one excitation with the smallest 
energy consists of removing a spin-down particle from Landau level and adding a spin-up 
particle to Landau level 1, implying a mean-field spin gap of huOf.. For any momentum Q, 
the imaginary part of the spin susceptibility Imx(Q,t^) consists of delta functions at the 
frequencies uj equal to multiples of the cyclotron energy uoc. The mean-field spin susceptibility 
is represented as the thick vertical lines in Fig. 0, in which, for Q = 0.6/ao, the height of 
each line is proportional to the relative strength of its delta function. For smaller Q, the 
delta function ai uj = ujc dominates the others to an even greater extent. 

Our final result for the spin susceptibility of the fractional-statistics gas differs substan- 
tially from the mean-field value. Fig. Q displays the imaginary part of the spin susceptibility. 



calculated in Section as a function of the frequency to for various values of the momentum 



Q. For any momentum, the lowest energy excitation has an energy A of 0.32 ^cUc- Actually, 
the calculation has rather large error bars, as discussed in section |V|, and we only assert that 
the spin gap is less than 0.4 hujc- A numerical approach is needed to determine the spin gap 
more precisely. In any case, the spin gap, if it exists, is substantially less than the mean-field 
value of huJc- The susceptibility spectrum has broadened considerably, though half of the 
spectral weight lies within 0.01 ujc of the spin gap A. At larger momenta, a definite peak 
appears at the frequency u = 1.2 Uc- 

Interactions with phonons, the density oscillations of the system, broaden the spin sus- 
ceptibility spectrum into the form displayed in Fig. |^. Consider the effect of the phonons 
upon the single-particle density of states of the unoccupied Landau levels n > 1. As illus- 
trated in Fig. 0(a), a particle in Landau level 2 may emit a phonon and decay into Landau 
level 1, resulting in the broadening of Landau level 2. Also, assuming Landau level 1 is 
broadened, as in Fig. 0(b), a particle in that Landau level may make an intralevel transition 
to a lower energy state within the Landau level and emit a phonon, in turn broadening 
Landau level 1 self-consistently. This broadening extends the particle spectrum to lower en- 
ergies. Similarly, the intralevel transitions within Landau level broaden the hole spectrum, 
so that it reaches higher energies. Thus, the energy gap separating the particle and hole 
states decreases by an amount of order one, since the theory contains no small parameters. 
The phonon interaction between the particle and the hole lowers the spin gap further. 

Conceptually, a flipped spin modifies its local environment to diminish the overall exci- 
tation energy, the spin gap. In other words, the lowest energy state |q> with spin one and 
momentum q may be approximated by density waves in addition to spin waves, 

$> . (L13) 

Here |<l>> is the ground state, is the Fourier transform of the spin-raising operator cr~^{r) 
defined in Eq. ( |1.6|) , 

a+ = y c/ra+(r)e~*'i-'" , (L14) 



k 
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and the density operator pq is defined in Eq. (|1.21|) . The coefficients must be determined, 
and we have left out terms involving more than one phonon. In real space, the excited state 
|q> may be written. 



i 

where h{r) is the Fourier transform of h]f_. 



;i.i5) 



r 

Mr)=i^7^/^ke^"-, (1.16) 
CTj^ is the spin-raising operator on particle i, and Z has the value. 

The other particles are repelled from the flipped spin present in a higher Landau level. At 
the level of a Hartree approximation, the reduced density near the flipped spin leads to a 
reduced mean gauge field and a reduced cyclotron energy huc, and thus a reduced excitation 
energy. In conclusion, the single-mode approximation of the spin-one excitation made by 
Girvin et. al, [g 

|qsMA>= fx+ |$> , (1.18) 

describes the lowest lying spin-one excitation poorly. 

Let us define certain conventions used throughout this paper. We shall often use units 
in which the cyclotron energy huc, defined in Eq. ( p..l2|) , and the magnetic length oq. 



ao = ^=[2nil-u)p]-y' (1.19) 
both equal one. Let us define the density operator p(r), 

N 

pir) = J25iv-r,), (1.20) 

1=1 



along with its Fourier transform pq. 



N 

P^ = ll^-'''-'' ■ (1-21) 
1=1 



The mean-field current operator j(r), 

^ 1 r e- 1 
jW=E^{P^ + -A.,5(r-r.)} , (1.22) 

has a Fourier transform jq, 

J, = gi{P. + ;A.,e— }, (1.23) 



and a transverse component, 



jj=jq-(zxq), (1.24) 



where A represents the constant magnetic field of Eq. ( |1.10 ) 



II. EFFECTIVE INTERACTION 

Interactions with phonons, the density oscillations of the fractional-statistics system, 
broaden the spin-susceptibility spectrum and lower the spin gap. To analyze the system, we 
shall replace the Hamiltonian Tianyon? Eq- ( |1 ■ 1| ) , with an effective Hamiltonian consisting of 
particles in a magnetic field coupled to phonons. The phonon system accurately describes 
the spin-one excitation spectrum of the fractional-statistics gas. Perturbative techniques 
will then determine the spin susceptibility and spin gap. 

The effective Hamiltonian Tics modelling the fractional statistics gas has the form 

''^ 1 e ^ 1 ct 

^eff = 5^ 75— Pi + - A, + ^ -hug (aj^aq + aqajj) + ^ — i= jj (a\ + a^) . (2.1) 

i=l C q Z ^ ' q ^/ZU!q 

The operators and Oq create and destroy a phonon of momentum and energy huq. 
The particles lie in the constant magnetic field A given by Eq. ( |1.10| ), arising from the mean- 



field Hamiltonian Ti^, Eq. ( |1.9| ). The mean-field transverse-current operator jq, defined in 
Eq. ( [1.24| ), couples a particle to a phonon with coupling strenth aq. As in the mean-field 
Hamiltonian Tig, the particle density corresponds to filling the lowest Landau level due to 
the magnetic field A with particles of both spins. The phonon dispersion Uq and coupling 
aq, along with the justification of the effective Hamiltonian Ti^g, are given below. 



Let us first approximate the part of tlie Hamiltonian left out of the mean-field treatment 
as a two-body interaction neglecting the three-body terms which have a limited effect 
upon the spin gap. The perturbation Hamiltonian Tii, defined by subtracting the mean- field 
Hamiltonian, Eq. ( |1.9| ), from the full fractional-statistics Hamiltonian, Eq. ( p..l[ ), in the limit 
of large a, may be written, 



T^l — Tianyon — T^Q 

The potential V(q) is defined to be, 

V(q) = 



(2.2) 



^1 . I he 
'l-jy) — — (z X q) 



iL"^ eq 



(2.3) 



and the density operator pq and the mean-field current operator jq are defined in Eqs. (|1.21|) 
and ( |1.23| ). The second term of Hi, the three-body interaction, is largest for p equal to q, 
where the density operator pp_^ reduces to the total particle number N. Removing the terms 
with p not equal to q leaves the two-body interaction Hamiltonian, 



^2-body ^ ^ ^ ^(q) ■ j^Pq - ^ r E |V(q)lVqP^ 



1 



V°(q) 



Pq 
Jq 



(2.4) 



where, for q in the x direction, the bare interaction V° has the form. 



1 \) iq 



V 



(2.5) 




-iq 

In the preceding equation and the remainder of this paper, the cyclotron energy hujc, 
Eq. ( |1.12D , and the magnetic length oq, Eq. ( |1.19D , are both set equal to one. 
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Various higher order processes modify the effective interaction between any two particles, 
in a fashion similar to the screening of the Coulomb interaction in an electron gas. Dia- 
grammatically, the effective interaction is the sum of a bare interaction and a polarization 
bubble attached to two bare interactions, 

V(q, u) = V°(q) + V°(q)P(q, a7)V°(q) . (2.6) 

The full polarization bubble V is the correlation function 



J in Jo 



X <<l>| 



^r{p(r,t),p(0,0)} T{p(r,t),j-(0,0)} T{p(r, t), jJ'(0, 0)} ^ 
r{j-(r,t),p(0,0)} T{j-(r,t),j-(0,0)} T{f{r,t), jy{0,0)} 
\T{jy{r,t),p{0,0)} T{jy{r,t),f {0,0)} T{j^(r, t), j^(0, 0)}^ 



|$> . 



(2.7) 



The density operator p(r) and the mean-field current operator j(r) are given by Eqs. ( |1 . 20|) 



and (|1.22|) , the time dependence of the operators follows the definition of the Heisenberg 
representation, Eq. (|1.4]) , and the operators are time ordered. Eq. (|2.6|) for the effective 
interaction may be expressed in terms of the proper polarization V^, the Feynman diagrams 
in the polarization which remain connected after removing any single interaction line, in the 
form, 

V(q, Lu) = V°(q) + V°(q)P^(q, ^)V(q, u) , (2.8) 

as displayed in Fig. H(a). The bare interaction V*^, Eq. (|2.5|) , is denoted by a thin wavy line, 
the effective interaction V is denoted by the thick wavy line, and the proper polarization 
is denoted by the black square. 

Let us calculate the effective interaction using the Random Phase Approximation. In the 
Random Phase Approximation (RPA), the proper polarization in Eq. ( |2.6| ) is estimated 
by the bare polarization of the mean-field Hamiltonian TYq, Eq. ( |1.9D . As displayed 
diagrammatically in Fig. §(b). 



V^^^(q,a;) = V°(q) + V°(q)P°(q,a;)V^P^(q,u;) 



(2.9) 



where the straight hnes represent single-particle propagators of the Hamiltonian TYq- The 
bare polarization D° is found by replacing |$> with |$o>5 the ground state of TYq, in 
Eq. ( ^.71 ) and using the Hamiltonian T-Cq in the Heisenberg representation, Eq. ( |1.4| ). The 
calculation of the bare polarization of a spinless system in and is modified for spin | 
semions to have the value, 



P°(q,u;) = 2 



27r 



(2.10) 



where 



-{n-by 



(2.11) 



^1 (^^-l)!K - {n-ivY] 
and 6 represents The factor of two in Eq. ( p.lO| ) arises from the spin degrees of freedom. 
Substituting Eqs. ( p.5| ) and ( |2.1CI| ) into Eq. ( p^Ql ) results in the RPA interaction. 



RPA, 



1 27r 



2 L2g2£) 



I + S2 ig(l + Si)^ 


-^g(l + Si) q^Jlo 



(2.12) 



where 



D 



:i + si)2-so(i + s2 



(2.13) 



For fixed q, the RPA interaction V^^^{q,uj) contains poles at discrete values of u, with 
a spacing of approximately one between poles. Fig. |^ graphs the imaginary part of the 
current-current component of V^^^ for three values of q. The solid line in each plot has 
been broadened with a small value of rj for visibility. However, that discrete Landau level 
structure is not physical, and would disappear in a more accurate calculation of the proper 
polarization V^. The effective interaction at finite q will be broadened over a range of 
energies, resembling the dashed line in Figs. |^(b) and|^(c). 
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Nevertheless, we will approximate the dynamic effective interaction as a sharp mode for 
each momentum q in the current-current channel, disregarding both the broadening of the 
effective interaction and the effect of the density-density and density-current channels. For 
g < 2, the lowest energy singularity in the current-current RPA interaction dominates the 
other singularities in the three interaction channels, while for q > 2 the dynamic interaction 
is too weak to have a significant effect. The effective interaction then has the form. 



V'="(q,cu) = V"(q) + 














\ 



(2.14) 



The excitation, a phonon, has an energy liLOq and a transverse-current coupling a^. 

Let us determine the coefficients and lo^ appearing in Eq. ( p.l4| ) for the effective 
interaction V*^^. We will define a quantity equal to the area under the curves in Fig. ^ 

1 



ImVi^;^(g,a; 



dijj . 



(2.15) 



Then ujq is defined as the weighted average value of a;. 



1 



vr 



and the coupling strength a„ obeys 



[2.16) 



\aq\ =2AqUq. 



[2.17) 



Substituting Eq. (|2.5|) into Eq. (pl8|) demonstrates that 



2 PP ' 



which, along with the f-sum rule [|r^ , reduces Eqs. ( P^.16D and ( ^.17] ) to 



(2.18) 



TT 



an 



L2 • 



(2.19) 



Numerical calculations confirm this constant value of la^P, and fit oOq to the formula. 
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(formula) 
9 



VsqY + (gV2)^ , (2.20) 



with a sound velocity Vs of 1. Both tUq and 

(formula) 

graphed in Fig. 0. 

Actually, the Random Phase Approximation does not determine the speed of sound Vg 
sufficiently accurately. Together, Hartree-Fock and correlation effects [^^3 lower the sound 
speed of the spinless fractional-statistics gas by approximately 9%. The sound speed of the 
spin i fractional-statistics gas would presumably shift by a comparable amount. Inserting 
Vs = 0.9 into Eq. ( p.20|) has a sizable effect upon the resulting spin gap, which we included 



in the error bars of the spin gap calculation in Section 0. 

The static interaction V° in the effective interaction V^^, Eq. ( p.l4|) , does not affect the 



spin gap. When calculating the spin susceptibility diagrammatically, as in Fig. 0, each of the 
phonon interaction lines is short-ranged. Dampening the interaction Tianyon — '^^o for large 



particle separations does not affect those Feynman diagrams. By Kohn's theorem [|Tl|], the 
remaining short-ranged interactions between particles in a magnetic field do not change the 
excitation energy at small momenta from the value without an interaction, the cyclotron 
energy huc- Thus, the static interaction V° does not change the Q ^0 spin gap and may 
be omitted. Beran and Laughlin demonstrated that with a cutoff parameter a in the 
fractional-statistics Hamiltonian 'H^^-^y^^-^-, Eq. corresponding to a reasonable doping of 

a lattice model, the static interaction V° does not create a roton to lower the spin gap at 
finite Q. Thus, the fractional-statistics system may be described as the mean-field Hamil- 
tonian Tip, Eq. ( |1.9| ), and a dynamic transverse-current interaction, the second term in 



Eq. ( p.l4| ), all of which reduces to the effective Hamiltonian IH-^^, Eq. (|2.1|). 

In summary, we have obtained a simpler system 7^^^, Eq. ( |2.1| ), which possesses the same 
spin gap as the fractional-statistics gas H^^^^^, Eq. (p..lD. The effective Hamiltonian consists 
of particles in a magnetic field coupled transversely to phonons. The phonon dispersion Uq 
is given by Eq. ( p. 20 ), and the particle-phonon coupling aq is given by Eq. (|2.19|) . In the 



remainder of this paper, we will investgate the effective Hamiltonian Ti.^^ to find the spin 
gap. 
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III. SINGLE-PARTICLE DENSITY OF STATES 



Let us compute the single-particle density of states of the fractional-statistics gas. The 
density of states allows a crude estimate of the spin gap, since interactions between a particle 
and a hole present in a spin-one excitation lower the spin gap further. However, the single- 
particle computation clearly illustrates the mechanisms of broadening and gap lowering 
present in the spin susceptibility spectrum. Also, the formal methods used to compute the 
spin susceptibility in section |V| resemble the techniques used to compute the density of 
states. We will start with the effective Hamiltonian 'H^^, Eq. (|2.1| ), consisting of particles 
in a magnetic field coupled to phonons. As displayed in Fig. ^, a particle in any Landau 
level may decay into either a lower Landau level or a lower energy state in the same Landau 
level, emitting a phonon. In this section, we will find the extent of the spectrum broadening 
induced by that decay process. 

The density of states depends upon the single-particle Green's function of the fractional 
statistics gas, 

G(ritio-i|r2t2a2) = -z < $ I T{^,,(ri,ti) ^>l^{r2M)} I $ > , (3.1) 

where |$> is the ground state. The operators \E'J.(r,t) and \E'o-(i',^) create and destroy a 
spin (J particle at position r and time t in the Heisenberg representation, Eq. (|1.4| ). The 
particle may emit and later absorb a virtual phonon, giving the Green's function a complex 
self-energy which broadens the spectrum and lowers the spin gap. As illustrated in Fig. ^ we 
will solve Dyson's equation for the exchange of a single phonon in the effective Hamiltonian 
Ti^g, Eq. (|2TT|), to find the spin gap. 

Let us first find the Green's function corresponding to the bare Hamiltonian Eq. ( p.. 91 ), 
in the absence of phonons. A constant magnetic field in the symmetric gauge possesses the 
Landau level orbitals ipnk [H], 



k 

|2/ 



^^,(^) = ^^^_(i^_2A) (1^*-2|-1 e-l^l^/S (3.2) 

^ ^' V27r2"+fcn!A;! V2 dz* ) \2 dz J ' ^ ^ 

with Landau level energies, 
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eS=ef =n+l. (3.3) 



The complex number z = x + iy refers to the position, and the natural units hujc, Eq. ( 1.12 ) 



and ao, Eq. ( 1.19| ), are set equal to one. As drawn in Fig. |I], Landau level is occupied with 



particles of both spins, while the higher Landau levels are empty. The bare Green's function 
may be written in the Landau level basis as, 

G°(riti(Ti|r2t2a2) = S<r,a,^J E E ^' I C^^) \nk> Vn'k'{v^Wnkij2)e'''^''-'''^ dE , 

nk n'k' 

(3.4) 

where 

<n'k'\ G\E) \nk>= <5„„,4,,^^-^-^^_ , (3.5) 

E - [n + 1/2) + ir]n 

and the infinitesimal ?7„ is positive for levels n > 1 and negative for level n = 0. 

The Green's function of the effective Hamiltonian Ti^g, Eq. ( p.l|) , may be found by solving 
Dyson's equation for the exchange of a single phonon. This process is illustrated in Fig. |^, 
where the thin straight line is the bare Green's function G*', the thick straight line is the 
full Green's function G, and the thick wavy line is the phonon propagator. We included the 
diagrams in which, besides the phonons, only a single particle or a single hole is present at 
any point in time. For instance, the exchange graph in Fig. |] where n > 1 and m = has 



been neglected. This approximation is similar to the Tamm-Dancoff approximation in 
nuclear physics. We have also left out diagrams with crossed phonon interaction lines. The 
relevant particle quantum states \n k>, in addition to the phonons, consist of Landau level 
filled with particles of both spins, with an extra particle in the orbital (fnk for n > 1, or a 
hole in the orbital ipok for n = 0. 

Dyson's equation involves the matrix elements of the transverse-current operator j^, 
Eq. ( |1.24| ), between the Landau level states \n k>. Summing over the intermediate states 
\mp> within a Landau level m and averaging over the direction 6 of the phonon momentum q 
combines those matrix elements into a simple expression. After performing those operations, 
we find, 

14 



1 r^^. 

271 Jo 



/ V <n' A;' I I m]9><mp I jTq I n /c> (i6' = 5„n'^fcA:'|A<nm(g)|% (3.6) 
Jo p ^ 

The effective matrix element M.nm{Q), which is independent of the index k, has the value 



Mnm{q) = \/^&^""""'^^' [nLl-^-\b) - e'^/^ , (3.7) 

where L denotes an associated Laguerre polynomial, 

L^ib) = yb-'^— (e-^6"+") , (3.8) 

and b stands for 

We may now write down Dyson's equation for the single-particle propagator arising from 
the exchange of a single phonon. The Green's function will be computed in the Landau level 
basis, 

G(ritiai|r2t2^2) = S,,,,^ /EE <^'^' I G{E) \ nk> ipn,u'{r^WnA^2)e'''^''-''UE . (3.9) 

nk n'k' 

Because of Eq. ( p.6|) , the full Green's function G{E) is diagonal in the basis of Landau level 
orbit als, 

<n'k'\ G{E) \nk>= Snn'Skk'Gn{E) . (3.10) 
The Landau level Green's function G„ has the standard form, 

E - {n + 1/2) - T.n{E) + ir}n 

with a self-energy, 

E„(E) = ,L^rf^/^i:^l=M!KI!G„(E + .). (3.12) 



The coupling la^P is given by Eq. ( p.l9| ), and the phonon energy ujq is approximated by 



Eq. ( p.20| ). When n is an unoccupied Landau level, n > 1, the index m is summed over the 
levels m > 1, and when = 0, only m = is present in the sum. 

The analytic structure of the Landau level Green's function Gm{E) allows a simplification 
of the self-energy expression, Eq. ( |3.12|) . When Landau level m is unoccupied, each of the 
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poles of Gm{E + uo) lies in the lower half of the complex plane of cu, while the phonon 

propagator has a pole on each side of the real axis, as indicated in Fig. ^. The uj integral 

may be evaluated by closing the contour upwards, so that for n > 1, 

r da °° \n P 

^n{E) = L' \MUq)r-^GUE - u,) . (3.13) 

l^^J m=l 

Landau level has a similar self-energy expression. 

Even though the resulting density of states is physically incorrect, it provides a reasonable 
estimate of the spin gap. We solved Eqs. (PTTTD , ( PTTB[ ), and ( PTTID for Gn{E). In Fig. [T^, 



we display the density of states Dn{E) within the lowest four Landau levels, 

1 



Dn{E) = p 



ImGniE) 



(3.15) 



Our calculation left out the effect of the static interaction, Vq in Eq. ( 2.14| ), which changes 



the density of states greatly. The static interaction causes a gap between the occupied and 
unoccupied states which diverges logarithmically with the size of the sample. However, due 
to Kohn's theorem, as argued in Section 0, the static interaction does not affect the spin 
gap, the quantity of physical interest. Thus, except for leaving out the particle-hole and 
multiple-phonon interactions, the resulting spin gap remains valid. 

The resulting broadening of the spectrum significantly reduces the gap energy relative to 
the mean-field gap displayed in Fig. The broadened unoccupied states in Fig. [1^ extend 
down to a minimum particle energy of 

^min.p = (ef ^ - 0.48) hcoc = 1.02 fiujc , (3.16) 

and the occupied states reach a maximum hole energy of 

^max,h = (eP + 0.05) huj, = 0.55 hcu, . (3.17) 



If the particle-hole interactions are neglected, the lowest lying spin-one state then has an 
energy, 
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A- P-^ = ^^i„,p - E^ax,h = 0A7huj, . (3.18) 

Both the dynamic interaction between the particle and the hole and multiple phonon inter- 
actions will lower the spin gap further. 

Finally, let us examine the overall distribution of the single-particle spectrum. The 
density of states of Landau level 1 diverges as the energy approaches -Emin.p from above, so 
that roughly 60% of that level's states are within 0.01 of -Emin,p; while the tail extending 
to higher energies contains the rest of the states. Similarly, 90% of Landau level lies 
very close to -Emax.hj with the tail containing the remaining spectral weight extending to 
lower energies. Landau level 2 contains a strong peak at an energy somewhat below e^2 \ 
addition to some broadened structure at both higher and lower energies. Landau level 3 is 
more broadened than the lower Landau levels, though a peak remains visible. Despite the 
broadening, the average energy of each Landau level, 

1 r°° 

<En>=- EDn{E)dE, (3.19) 

P J — oo 

maintains the unperturbed value ej^\ 

IV. SPIN SUSCEPTIBILITY AND GAP 

Let us compute the spin susceptibility of the fractional-statistics gas, and consequently 
its gap. We shall find the retarded spin-one particle-hole propagator, 

J^(rir4ti I rarats) = -^^(^1-^2) <$ | ^I(r4,ti) ^T(ri, ^1)^1(^2, t2)^i(r3, ^2) | $> • (4.1) 

The particle-hole pair produced in the intermediate state may emit and then subsequently 
absorb a phonon. The phonons arise in the effective Hamiltonian Ti^g, Eq. ( p.l| ), which as 
argued in Section |I| has the same spin gap as the fractional-statistics gas. The single-particle 
self-energy, S„(i?) in Eq. ( |3.13| ), accounts for the case when the particle both emits and 
absorbs a phonon, and the single-hole self-energy, So(-E') in Eq. ( ^.14| ), accounts for the case 
when the hole both emits and absorbs a phonon. However, the full particle-hole propagator 
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T is needed to include the sizable interaction between the particle and the hole, where one 
emits and the other absorbs a phonon. Consequently, the spin one spectrum exhibits more 
broadening and gap lowering than does a convolution of the single-particle Green's functions 
computed earlier. The spin susceptibility and spin gap are readily expressed in terms of the 
pair propagator. 



As in the section |TT[ we have only included processes in which a single particle-hole pair 
propagates forwards in time, emitting and absorbing phonons. The Random Phase Approx- 
imation (RPA) Feynman graphs, which for the spin particle-hole propagator exhibit the 
corrections to the Hartree-Fock ground state most strongly, as a reversed time-ordering [0 , 



do not appear in the spin 1 propagator. Thus, the Tamm-Dancoff approximation [12| should 



be reliable. Graphs containing crossed phonon interaction lines have again been excluded. 



A magnetoexciton basis |jT3[ , describing a particle and a hole in a magnetic field in terms 
of eigenstates of both momentum and the bare Hamiltonian Ti^, proves to be convenient in 
studying the particle-hole propagator. The magnetoexciton particle- hole wavefunction has 
the form, 

where z = x + iy is a. position written as a complex number. The spin-up particle lies in 
Landau level n > 1, the spin-down hole lies in Landau level 0, and the pair has momentum 
Q, represented as a complex number Za = iQx — Qy The corresponding quantum state 
|n Q> is an energy eigenstate of the mean-field Hamiltonian Tig, Eq. ( |L9|) , with energy 

= ^ • (4-3) 

The bare particle-hole Green's function J-"^ of the Hamiltonian Ti^ may be expressed in the 
magnetoexciton basis. 



^°(rir4ti|r2r3t2) = :^ /E E Q' I -^°(^) I ^ Q> 



n Q n' Q' 

X Vn'a'(ri, r4) Ca(r2, rs) e*^^*^-*^) dE , (4.4) 
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where 



1 



(4.5) 



E — n + irj 

The phonon couphng appearing in the particle-hole propagator depends upon the 
transverse-current matrix element between two magnetoexcitons. The transverse-current 
operator j^, defined in Eq. (|1.24|) , may act on either the particle or the hole, and both 
couplings must be included. The coupling between the particle and a phonon has a matrix 
element, 



<m Q-q I jj^^^ I n Q> 



X 



-(Pi-PiO + A(ri) 



dvidY2 lim e ^'^'"^ 
■ eq C.fl(r'i,r2)^/'na(ri,r2) 



while the coupling between the hole and a phonon has a matrix element. 



(4.6) 



<m Q-q I jq I n Q> 



dri dr2 lim e 



X 



fP2-P2')-A(r2) 



eq ^m/3(ri.r2)^„a(ri,r;,) 



(4.7) 



The mean-field A(r) is specified by Eq. ( |1.10|) , eq is the transverse direction z x q, and 
Zjs = i{Qx — (lx) — {Qy~<ly)- Siucc the hole has the opposite charge, its interaction with both 
the phonons and the magnetic field has the opposite sign. The resulting matrix elements 
have the values, 



and 



X 




V2, 



(4.8) 



(4.9) 



where = iq^—qy, b denotes and the Laguerre polynomials L are defined in Eq. ( p. 81) . 
The transverse-current matrix element is the sum of the particle and the hole components. 
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<m Q-q I jq I n Q> = <m Q-q \ jq^^UnQ> + <m Q-q | jq | n Q> 



X l^e^im^'^xQ)-^ [nL;;-'"-i(6) 

- ^nme^^'/'^^'^'''^^" • (4.10) 

Since small momenta dominate the interaction, we may neglect the phases proportional to 
q X Q in Eq. ( |4.10| ), approximating the current matrix element <m Q — q | Jq | ^ Q> by a 
term A1^^(q) independent of Q, 

MU^) = (^j [nLl—\b) - bLl--+\b) - e-'/' , (4.11) 



which ioT n ^ m has the same magnitude as A^nm(q), defined in Eq. ( p.6| ). 

The Hamiltonian Ti^g, Eq. ( p.l| ), restricted to a single particle- hole pair being present, 
may be treated exactly as a discrete quantum state coupled to phonons, allowing the dia- 
grammatic computation of the particle- hole propagator JF, Eq. (|4.1|), and the spin gap. Since 
momentum is conserved, consider states in which the system has total momentum Q. The 
particle's Landau level n and the set of phonons present specify a basis for the single-pair 
configurations. A system which may occupy one of a discrete set of quantum states \n>, 
where n is a positive integer, interacting with phonons has an isomorphic Hilbert space. If 
the corresponding matrix elements are the same, the spectrum of the discrete system agrees 
with the spectrum of Ti.^^ constrained to the restricted Hilbert space. Thus, Ti.^^ with a 
single particle-hole pair may be transformed into the discrete system. 



OO 



n=l q 2 



+ E E E ^ M'^M^l^n K - «-q) • (4-12) 

n=lm=l q V^^? 

The operators \l'Jj and \1'„ create and destroy a particle in the discrete state \n>. 

We found the particle-hole propagator by solving Dyson's equation for the exchange of a 
single phonon in T^^j^^^t^, as illustrated in Fig. |TT[ In effect, this method includes the particle 
exchange, hole exchange, and particle-hole ladder Feynman diagrams. The bare particle-hole 
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propagator JF°, defined in Eq. ( [4. 41) , is represented by a thin dashed hne, while the phonon 
propagator is denoted by a thick wavy hne. The fuU pair propagator T ^ represented by a 
thick dashed hne, may by examined in the magnetoexciton basis, 

^(rir4ti I rsrats) = :^ / E E Q' I •^(^) I ^ Q> 

n Q n' Q' 

X V^„,,,(ri,r4XJr2,r3)e*^(*^-*^)t/E . (4.13) 

Due to the phases of A^^^(q) in Eq. (^4.11|) , the integral over the direction of the phonon 
momentum q cancels in Dyson's equation whenever n ^ n', keeping the propagator diagonal. 

The resulting Dyson's equation of the particle-hole Green's function has a form similar 
to Eqs. (|3.1CI|) through ( |3.12| ) for the single-particle Green's function. The particle-hole 



propagator remains both diagonal and independent of momentum in the magnetoexciton 
basis, 

<n'Q!\ T{E) I n Q >= 5nn' 5qq' -^n(^) • (4. 14) 

The magnetoexciton propagator J-'n{E) for particle Landau level n obeys the equations, 

TniE) = ^ .p., . , (4.15) 

E -n~ 2ln[E) + IT] 



with a self-energy, 

.^2 r duj f ^ \M',,rn{q)?\aq? ^ , . \ (A^a\ 

^"(^)=^^ loo 2^ y (2^ S - - ^r^f "^^ ^ ■ ^ ^ 

Performing the integral over cj, using the analytic properties of T identical to the analytic 
properties of G leading to Eq. (|3.13| ), simplifies the self-energy to 

^n{E) = L' / 7^ E \MU<l)r-7^ ^m{E - UJ,) . (4.17) 



Eqs. (|iA5| ) and ( ^T7|) may be solved for J^„(E). 



Let us compute the spin susceptibility x(g, cj), as defined in Eq. (|1.3|), and thus find 
the spin gap. Setting ri equal to and r2 equal to in the particle-hole propagator 
JF, Eq. ( [4.13| ), and Fourier transforming yields the spin susceptibility. In terms of the 
magnetoexciton basis, 
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x(Q,w) = 7^ E l<0 I PQ I Q>|'-F„H , (4.18) 

where < | pq | n Q > is the density matrix element between the magnetoexciton |n Q > and 
the ground state, 

<0|pQ|nQ>= |rfri^L(ri,rOe-'i- = -^(^^j e-"\ (4.19) 

where b = \Q'^ and = iQx—Qy The resulting spin susceptibility is graphed in Fig. || and 
discussed in the introduction. This diagrammatic approach finds a spin gap A of 0.32 huJc- 
The next section will apply a different computational technique to estimate the error bars 
of the spin gap calculation. 



V. SECOND-ORDER PERTURBATION THEORY 

The spin gap may also be investigated with second-order perturbation theory. We can 
thereby estimate the uncertainty arising from leaving out the Feynman diagrams containing 
crossed phonon lines. Since the coupling constant aq is ^/7T in the natural units of the 
problem, processes involving multiple phonons have a noticable effect, lowering the spin 
gap further. In fact, the lower spin gap energy given by second-order perturbation theory 
appears to be more accurate that the value arising from summing diagrams, as discussed at 
the end of Section |VT[ However, second-order perturbation theory cannot determine the full 
spin susceptibility spectrum, but only the spin gap. 

Let us compute the spin gap, the lowest energy spin-one state, using second-order per- 
turbation theory. For a system with a ground state |0> and excited states |/> possessing 
energies E^^^^ according to a bare Hamiltonian Ho, an interaction Hi shifts the ground state 
energy to the value 

^ -^0 2^ (0) • y^-^) 

We will find the energy shift of the spin 1 particle-hole state |1 Q>, which serves as the 
ground state |0> in Eq. (|5.1| ). The magnetoexciton |1 Q>, defined by the wavefunction in 
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Eq. ( [4.2|) , has a spin-up particle in Landau level 1, a spin-down hole in Landau level 0, and 
total momentum Q. The relevant excited states |/> consist of the phonon of momentum q 
and the magnetoexciton |m Q — q>. 

Inserting these states and the interaction of Ti^^, Eq. ( |2.1| ), into Eq. ( |5.1| ) leads to the 
spin gap, 

\{2) _ M (o)\ 2 / "vL Fgl I ''t h i jq i ^ 

Replacing the current matrix element in Eq. ( ^.21) with M.'^i{q), as defined in Eq. ( [4. Ill) , 



and using the values in Eq. ( pl3|) for e^°^ yields the lowest energy of a spin-one excitation, 

A(^) = 1 - / A ^ E ^^^^^^ = 1 - 0.83 . (5.3) 
J (27r)2 2ug ^^m-l + Uq ^ ^ 

Thus, at the level of second-order perturbation theory, the phonon interaction reduces the 

spin gap to 0.17 huc, less than 20% of its mean-field value, compared to a spin gap of 0.32 hujc 

obtained by diagrammatic methods in Section 

Altering the sound speed Vg in the phonon dispersion, Eq. ( |2.20| ), changes the energy 

gap significantly, increasing the error bars on the spin gap calculation. For instance, a 10% 

decrease in the sound speed to Vs = 0.9 reduces the second-order perturbation theory energy 

gap calculated in Eq. (|5.3|) to 

^(2) ..=0.9 _ ^ _ Q_gg _ (5 4) 

Since such a sound speed is plausible 0, we cannot set a lower bound on the size of the spin 
gap, or even prove that a spin gap exists. Still, the spin gap definitely lies in the range 

A = {0.2±0.2)hujc, (5.5) 

much less than the mean- field value of hoOr- 



VI. BOSON TEST CASE 



We will compute the Green's function of the non-interacting spinless bose gas in order 
to test the computational procedure used to find the spin gap of the spin ^ fractional- 
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statistics gas. The same gauge transformation changing particles with fractional statistics 
into fermions obeying Hanyon, Eq- (P^), will also transform a gas of non- interacting bosons, 

N IP 12 

Wb.. = g^. (6.1) 

into Tianyon with a statistical factor z/ = 0, and with the cutoff parameter a—^oo. Since the 
bose energy eigenstates are known exactly, the bose Green's function may also be computed 
exactly. Agreement between the exact and approximate bose Green's functions demonstrates 
the validity of the semion calculation, because bose statistics involve a larger deviation from 
non-interacting fermions than do fractional statistics. 

We shall compute the part of the zero-distance bose Green's function, 

/oo 
^iEt ^ ^ I 2-{^(o, t) (0, 0)} I $ > rft , (6.2) 
-oo 

arising from the occupied states, where |$> is the ground state. In the fermionic repre- 
sentation, \l/'''(r) and ^^(r) are local fermionic operators creating or destroying a fermion 
at position r. In the bosonic representation, however, those operators also multiply the 
wavefunction by a phase depending upon the positions of all the other particles, effectively 
creating or destroying a vortex. Thus, the Green's function defined in Eq. ( |6.2| ) differs sub- 
stantially from the traditional bose propagator. In terms of the excited states |/ > with 
energies e/, the Green's function has the value, 

|</ I ^t(o) I $> ' |</ I $(0) I $>|' 



GiE) = y: 



(6.3) 



E — ei + if] E + ei — if] 

The imaginary part of the Green's function corresponding to the occupied states is then, 

- lmG{E) = V |</ I ^(0) \<^>f6{E + q) , (6.4) 

valid for E < 0, which are the only energies we shall consider. 

The Green's function calculated with the approximation techniques used to find the 
spin gap of the fractional-statistics gas agrees fairly well with the exact boson spectrum, 
as displayed in Fig. |12[ We are using units in which oq, Eq. (|1.19| ), and Uc, Eq. ( |1.12D , 
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defined with z/ = 0, are both set equal to one. The exact bose Green's function, the solid 



line in Fig. |T2|, was found by substituting the exact bose wavefunctions and energies into 
Eq. ( p. 41 ). The approximate bose Green's function, the dashed line in Fig. |12|, was found 
by first determining the bose effective Hamiltonian in the form of Ti^g, Eq. ( |2.1| ). We 



then considered, in addition to the single-phonon exchange processes considered in Section 



nil , processes involving arbitrary numbers of crossed phonon interaction lines, such as the 



diagram in Fig. |T^. The effect of such crossed phonon graphs in the fractional-statistics gas 
is discussed at the end of this section. The exact and approximate Green's functions possess 
a similar overall structure. In addition, the highest energy hole states in the two calculations 
differ in energy by only 0.06 TiuJc-, demonstrating the reliability of the calculational technique 
used to compute the spin gap of the fractional-statistics gas. 

Let us compute the bose Green's function exactly. The system consists of bosons 
confined to a region of size L^, with an average density, 

N I , , 

The gauge transformation from H-^^^^, Eq. (|6.1|), to Ti^nyon) Eq. multiplies the bose 

wavefunction by a phase depending upon the direction of the separation vector between 
each pair of particles, but leaves the excitation energies unchanged. In the fermionic repre- 
sentation, the N particle ground state |$> has the wavefunction, 

i^>=:^nf^l, (6-6) 



where z = x + iy expresses a position as a complex number. The A^— 1 particle excited states 
|ki, . . . , kAr_i > consist of plane waves, with energies 

N-l 

i=i 

and wavefunctions 

N-l / _ \ (A^-l)! 

|ki,...,k^_i>= n W ^ e'^-'-^w+-+'^^'-^('v-i)] (6.8) 

i<j W^i T 
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summed over permutations r. The normalization A^, 

AT 



(6.9) 



L^-y(iV-l)!ni!---nM! ' 
depends upon the momenta kj with an occupancy Uj greater than one. 

The bose Green's function may then be calculated from the energy eigenstates. The 
matrix element in Eq. ( p.4| ) coupling the ground state |$>, Eq. ( p.6| ), to the excited state 
|ki, . . . , kiv-i>, Eq. has the magnitude, 



|<ki,...,kiv-i I ^(0) I $>| 



(6.10) 



Substituting these matrix elements and the excitation energies of Eq. (|6.7] ) into Eq. ( |6.4| 
leads to, 



^ImG(i.)^^(2.)^(--) 



E 



N-l 



{ki,...,k]v_i} 



"'I "'A'-l * -—1 



1 



1 / 

5|E + 5]-|k,f I dki---dk^_i. 



(6.11) 



-'1^1 I^N-l \ j=l 

Changing the sum into an integral generates a combinatoric factor cancelling out a similar 
term in the normalization A/". To avoid an infared divergence, we will cut off the momentum 
at a value fco, requiring a normalization of ImG'(-E') consistent with the density p. 



P 



dki ■ ■ ■ dk^^i 



(6.12) 



so that 



k - ^ 



With the cutoff, one can take the Fourier transform of lmG{E), 



(6.13) 



e''^'-lmG{E)dE ^ 

-oo 71 



\Jk>ko k^ 



N-l 



oo 1 > 



oo 1 N 

— de 



N-l 



^ - it{^ - 1) - it In ( -it^kl 



2^0 



(6.14) 
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Transforming back to E yields the exact hole Green's function, 
1 n r°° 

T™ /^exact / T7'\ r I 



1 j^^exact^^) = r e(V4)t[ln|t|-l]g-(V8)l*l ^-^{E-e•^^^^^)t ^ ^g^^g) 

TT 27r J-oo 



where 



er^* = -^ln(i?)-i[ln(2)-7] , (6.16) 
with = vri?^ and 7 denoting Euler's constant. The exact Green's function is graphed as 



the solid line in Fig. [12 



Let us now compute the hole propagator of the bosons with a technique similar to that 



used in Section 111 for semions, in order to test our approximation method in comparison 



with the exact result G^^^'^^lE), Eq. (|6.15|) . Some changes are necessary in the bose case. In 



particular, the collective mode of the bose gas has a quadratic dispersion, 

u;, = . (6.17) 

Even though the RPA calculation of Section ^ for a bose gas obtains a linear dispersion, we 
will use the quadratic dispersion, which would arise if all Feynman diagrams were included. 
However, the effective interaction V between two bosons in the fermionic representation will 
still be approximated by Eq. ( p.6|) , which was used in Section || to compute the effective 
interaction of the fractional-statistics gas. In that formula, the bare potential V° for bosons 
is twice as large as the semion case value, Eq. (|2.5| ). We will again approximate the effective 
interaction as the bare potential V° added to a sharp mode at each momentum q in the 
transverse current-transverse current channel, as in Eq. (p.l4|) . The polarization bubble V 



of the bose gas, defined in Eq. (|2.7|), satisfies the sum-rule argument in Section ||, resulting 
in a coupling strength of 

,9 27r 

\a,\' = — . (6.18) 



The static interaction V° moves each Landau level n from the bare energy e^^-*, Eq. ( |3.3| ), to 
its Hartree-Fock energy e^^, which for Landau level has the value 
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6^ = -^ln(i?) + ^[ln(2)-7] (6.19) 

for a sample of size L'^ = 'kE?. In summary, the resulting system is described by the effective 
Hamiltonian TieS, Eq- (|2.1|), with the phonon energy Uq given by Eq. ( |6.17| ) and the coupling 
Uq given by Eq. (|6.18| ). The energies corresponding to Landau level should be shifted by 
the amount e[^^ — eg'^^ 

Because the phonons have a quadratic dispersion, processes with crossed phonon lines, 
such as those in Fig. |13|, make a significant contribution to the bose spectrum. The solution 
to Dyson's equation for the exchange of a single phonon, Eqs. ( p.ll|) and (|3.14|) does not 



determine the propagator accurately. Instead, we will sum the Feynman diagrams with 
any number of crossed phonon lines. In the last paragraph of this section, we will discuss 
how such multi-phonon graphs affect the fractional-statistics gas. Since level is the only 
occupied Landau level, we shall only consider graphs where the hole lies in Landau level 
in the initial, final, and all intermediate states. 

To find the bose Green's function, a single Landau level interacting with phonons will 
be treated as a single quantum state interacting with phonons. The effective Hamiltonian 
Ti^g, Eq. (^T|) , representing the bose gas in the fermionic representation, Eq. ( |lT1) , will in 
turn be transformed into the discrete state Hamiltonian 

^single = + ttj^Oq) + ^ Ag^"^^ (a^ + ttLq) , (6.20) 

q ^ q 



similar to the Hamiltonian Ti.^^^^^^^^, Eq. ([4.12|) , utilized to find the particle- hole propagator 



of the fractional-statistics gas in Section |I^. The operators and create and destroy a 
fermion which may occupy a single quantum state, aj^ and Oq create and destroy a phonon, 
and the fermion-phonon coupling Xq has the value 



where Uq and aq have the bosonic values of Eqs. ( |6.17D and ( |6.18|) and is defined in 
Eq. ( p.6|) . The Green's function of Ti-^.^^^^, 
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/oo 
e^^(*2-*i) dE , (6.22) 
'OO 

where |$> is the vacuum state, may be converted into an approximate zero-distance bose 
Green's function defined in Eq. (|6.2|) , 

The differences between a discrete quantum state and a hole in any orbital of a filled Landau 
level account for the changes in sign and normalization. For processes involving the exchange 
of a single phonon, the self-energy expression for Ti^i^gi^ agrees with the single-hole self-energy 
formula Eq. ( p. 121 ) with n and m both set to Landau level 0. For the process involving two 



crossed phonon lines, as in Fig. |I3|, the spatial integral due to the Hamiltonian Ti^g , Eq. 
reduces to the sum of the matrix elements of transverse-current operator j^, Eq. ( L24 ) 
between the Landau level orbitals in Eq. (jsT 



^ <0 A;' |£q, I OpxOp |£q I Op'xOp' I jj I O/xOp" I jj | A;> 

= 4;c'e^('^'"'^)-'|^oo(g)n^oo(g )l' , (6-24) 



pp'p' 



Since small momentum phonons have the dominant effect, the phase e*^'^'^*'^'^ may be ne- 
glected. Then the two-phonon Feynman diagrams for both Ti^a and T^singie have the same 
self-energy expression, 

_ f°° du f dq duo' f dq' 



27r J (27r)2 J-oo In J (27r) 



|2 1^ |2 

X |^oo(g)riA^oo(g')r 



M|2 I'^pl \^q\ 



X Go{E + uj) Go{E + UJ + UJ') Go{E + u') . (6.25) 

This agreement may be generalized to Feynman diagrams with any number of crossed phonon 
lines. Thus, we may study H^i^gi^, to find the hole Green's function of the bose gas in the 
fermionic representation. 
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The Hamiltonian Tisingie, Eq. ( |6.2CI| ), whose Green's function is readily converted into the 
zero-distance Bose green's function, may be solved exactly by transforming the phonons, 
treated as harmonic oscillators, into the conjugate operators Xq^ and Pq at each momentum 
q. The phonon creation and annihilation operators a]^ and are the standard linear 
combinations of Xq and Pq, 

^P„ \ 



(6.26) 



and 



'KaXr 



/rria 



(6.27) 



where Kg and rriq satisfy 




(6.28) 



The Hamiltonian Hsingie may then be reexpressed as. 



"^single ^ ^ 

q L 

= E 



1 p2 ^ 

^Kgxl + 

2 ^ 'I 2mgj 



2K 



g1/ -q 
q V ^9 



(6.29) 



When the discrete quantum state is occupied, its interaction shifts the harmonic oscillators 
displacements and lowers the total energy by 



\2 



q 



(6.30) 



The vacuum state is then 



i$>=n 



e 2 



(6.31) 



with a length scale 



dg = {Kgnig)-^'^ 



(6.32) 



A relevant excited state contains the discrete state and nq phonons of each momentum q, 
with the shifted harmonic oscillator wavefunction 
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where the if„ are Hermite polynomials, and an energy. 



+ ^q\l 77 

y hnUJ. 



(6.33) 



(6.34) 



The excited states |nq>, Eq. ( p.33| ), and energies e„q, Eq. ( p.34|) , determine the imaginary 
part of the Green's function, defined in Eq. (|6.22|) , 



{"q} 



(6.35) 



The overlap between ^f^|$> and the excited state |nq> has a magnitude. 



where c„ stands for. 



q "q- 



(6.36) 



27T 



-e 2" 



(6.37) 



The energy shift AE in Eq. ( |6.3CI|) diverges due to phonons of small momentum. However, 
the state \E'"I'|$> has an expected energy of 



<^singlo>= AE + <r2q> Ug = 



(6.38) 



where <0> is shorthand for <$ | \& (9 | $>. Since on average the excited phonons cancel 
the energy shift, the energy of a state depends upon the fluctuations. 



= Z!(^q- <^q>)^g • 



(6.39) 



Substituting Eq. ( |6.36| ) and Eq. ( |6.39| ) into Eq. ( |6.35| ) and taking the Fourier transform leads 



to 



oo /I 



-ImG"^'°§^'=(E) 

TT 



{nq} q ^q- 



^—iiriq—Cq) Ulqt 



e " e 



exp \ ^ Cq [e'^^"^" - l+ituo, 
I q 

Jl/i) [(l+it) \n(l+it)-it\ 



(6.40) 
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where the imaginary part of the logarithm lies in the range, 

-| <Im[ln(l + zt)] < I . (6.41) 

After Fourier transforming back to energy dependence and shifting the t contour in the 
complex plane, the Green's function of TYgingie becomes, 

- - Im G'^-gi'=(E) = — e^/^-^ / e(^/^)*[i° 1*1-^] g-^^/^^l*! e^^* dt . (6.42) 

TT 27r J-oo 

The approximate zero-distance Green's function for the hole states, after applying Eq. ( |6.23D , 



- Im G™^(E) = Aei/4+(i^-.?^) r e(V4)i[in W8)l*l g-^C^-^D* , (6.43) 

TT 2tT J-oo 

is plotted in Fig. |1^ as a dashed line, and displays good agreement with the exact spectrum 

^cxact^^)^ Eq. ( |05D . 



The calculation of the spin gap in Section |V| left out the Feynman diagrams with crossed 
phonon lines. The method used in this section to deal with those graphs fails when Landau 
level excitations are involved. Although the higher Landau states affect the spin suscepti- 
bility, interactions involving only Landau levels and 1 cause a majority of the shift of the 
spin gap. If the Landau levels n > 2 are neglected, the preceding analysis of the fractional- 
statistics lowers the gap by the amount AE in Eq. (|6.3(]| ). Then, since AE is simply the 
second-order perturbation theory result found in Section 0, second-order perturbation the- 
ory gives the same spin gap as does including all the Feynman diagrams, for the Landau 
levels causing most of the gap's shift. Thus, the spin gap found by second-order perturbation 
theory is more accurate than the diagrammatic calculation in Section [rV| . 
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FIGURES 

FIG. 1. Illustration of the mean- field density of states of the spin ^ fractional-statistics gas. 
Landau level is occupied with particles of both spins. 

FIG. 2. Imaginary part of the spin susceptibility i^), defined in Eq. (|1.3D , at three momenta: 

(a) Q = 0.2/ao, (b) Q = 0.4/ao, and (c) Q = 0.6/ao, as calculated using Eq. (p^ ). The vertical 
lines (d) display the delta-function mean-field susceptibility at Q = 0.6/ao. 

FIG. 3. Illustration of decay processes causing broadening of the single-particle spectrum: (a) 
interlevel transitions and (b) intralevel transitions. Each decay process emits a phonon. 

FIG. 4. Diagrammatic representation of the effective interaction: (a) exact expression in terms 
of the proper polarization bubble T)^, as in Eq. ( |2.8| ), and (b) RPA approximate expression in 
terms of the bare polarization bubble V^, as in Eq. ( |2.9| ). 

FIG. 5. Imaginary part of Vjj, the transverse current-transverse current component of the 
RPA effective interaction given in Eq. (|2.12| ), at three momenta: (a) q = 0.5/ao, (b) q = I/oq, and 
(c) q = 4/ao. The solid lines are broadened with a small t] for visibility, while the dashed lines of 

(b) and (c) are broadened with a larger r] to resemble the full effective interaction. The collective 
mode energy is defined in Eq. ( |2.16| ) to be the average value of u>, and the coupling |agp defined 
in Eq. ( |2.17| ) is proportional to the integrated interaction strength at a given q. 

FIG. 6. Collective mode dispersion ujq. The solid line is the averaged RPA interaction defined 
by Eq. ( 2.16| ), and the dashed line is the approximation formula given in Eq. ( 2.20 ). 



FIG. 7. Typical Feynman diagram involving the static interaction, the V'^ of Eq. ( 2.14 ), which 
as discussed at the end of Section || does not affect the spin gap. 

FIG. 8. Illustration of Dyson's equation due to phonon exchange for the single-particle Green's 
function G, as expressed in Eqs. ( |3.9[ ) to ( |3.12| ). The labels on the particle lines refer to the Landau 
level orbitals given in Eq. (p.2|). 
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FIG. 9. Contour used to evaluate the u) integral in the self-energy expression E„(i?), Eq. ( ^.12 ) 
for Landau levels n > 1 and m > 1. 



FIG. 10. Density of states Dn{E), as defined in Eq. (3.15), for Landau levels to 3. The bare 
Landau level energies en^ given in Eq. (|3.3|) are indicated on the top of the graph. Neglecting the 



particle-hole interactions, the spin gap is the difference between the maximum hole energy Ej 



max.h 



and the minimum particle energy -Bmin.p- As discussed following Eq. ( 3.15 ), the calculation leaves 
out the divergent gap separating the occupied and unoccupied states. 

FIG. 11. Diagrammatic sum used to compute the spin one particle- hole propagator J^, defined 
in Eq. ( ^.1[ ). The particle- hole pair representation is transformed into a discrete system obeying 
^discrete' -^1- ( |4.12| ). The indices on the dashed lines refer to either the particle Landau level in the 
magnetoexciton wavefunctions given in Eq. ( [4. 21 ) or the states of Ti-^.^^^^^^ . 



FIG. 12. Green's function G, defined in Eq. (S.2), for the holes states of the non-interacting 
bose gas in the fermionic representation. The solid line denotes the exact spectrum G^^^'^^{E), 



given by Eq. ( |6.15 ), and the dashed line denotes the approximate spectrum G^^'^^°^{E), given by 



Eq. ( |6l3|) . The energy = on the graph corresponds to eg ^ , the Hartree-Fock energy of Landau 



Level given in Eq. ( |6.19| ). 

FIG. 13. Illustration of a process involving two crossed phonon interaction lines included in the 
bose gas calculation, as expressed in Eq. (|6.25| ). The hole propagators are labeled by the Landau 
level orbitals defined in Eq. ( |3.2D . 
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